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Abstract. Measurements of transverse-momentum (jpt) spectra of neutral pions in 
Au+Au and d+Au collisions at ^ii^ = 200 GeV and 62.4 GeV by the PHENIX 
' experiment at RHIC in comparison to p+p reference spectra at the same -^/snn are 

presented. In central Au+Au collisions at ^/snn = 200 GeV a factor 4—5 suppression 
for neutral pions and charged hadrons with pt > 5 GeV /c is found relative to the p+p 
reference scaled by the nuclear overlap function (Tab)- In contrast, such a suppression 
. of high-pj- particles is absent in d+Au collisions independent of the centrality of the 

C collision. To study the ^/snn dependence of the suppression Au+Au collisions at 

-i/snn = 200 GeV and 62.4 GeV are compared. 
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1. Introduction 
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The heavy-ion program at the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven 
National Laboratory is exploring the fundamental theory of strong interactions, QCD, 
under extreme conditions. A main goal is to create and study a deconfined and 
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thermalized state of strongly interacting matter, the quark-gluon plasma (QGP). Most 
particles created in central Au+Au collision at RHIC carry rather low transverse 
momenta (pr < 2 GeV/c). For the understanding of nucleus- nucleus (A+A) collisions, 
however, the small fraction of particles with high pj- originating from parton-parton 
interactions with high momentum transfer are of particular importance. In A+A 
collisions the hard parton-parton scatterings happen in the early stage of the collision, 
before a possible QGP has formed. The high-energetic scattered partons traverse the 
subsequently produced excited nuclear matter. As a consequence, particle production at 
high pt is sensitive to properties of the hot and dense matter created in A+A collisions. 

Here we study high-p^ neutral pions at mid-rapidity which are measured with the 
two central spectrometer arms of the PHENIX experiment. Each arm covers \rj\ < 0.35 
in pseudorapidity and A<f) = n/2 in azimuth. Neutral pions were measured by the 
PHENIX electromagnetic calorimeters (EMCal) via the tt° — > 77 decay. The EMCal 
consists of six lead-scintillator and two lead-glass sectors, each located at a radial 
distance of about 5 m to the interaction region pQ. The systematic uncertainty of 
the absolute energy scale of the EMCal is 1.5% in this measurement. 
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Figure 1. Nuclear modification factor Raa(pt) for neutral pions in central Au+Au 
collisions at */snn = 200 GeV compared to theoretical predictions as discussed in the 
text. 



2. High-p T Particle Yields in Au+Au collisions at v /5nn = 200 GeV 

Hadron production mechanisms in A+A are usually studied via their scaling behavior 
with respect to p+p collisions. Soft processes (j>t < 1 GeV/c) are expected to scale 
with the number of participating nucleons N part . Hard parton-parton interactions with 
small cross section, however, can be considered as an incoherent sequence of individual 
nucleon-nucleon collisions. In the absence of any medium effects the production of 
high-p T particles should be comparable to the production in p+p after scaling with 
a geometrical factor which reflects the increased number of scattering centers. It is 
customary to quantify the medium effects at high p T using the nuclear modification 
factor which is given by the ratio of the A+A to the p+p invariant yields [2] that are 
scaled by the nuclear overlap function (Tab)'- 

r U \ - d 2 Nf A /dydp T 

(T AB ) x d 2 a* p /dydp T 

The average nuclear overlap function (Tab) is determined solely from the geometry of 
the nuclei A and B. The average number of nucleon-nucleon collisions per A+B collision 
is given by (A coU ) = of^ x (T A b)- 

Raa{pt) measures the deviation of A+A from an incoherent superposition of NN 
collisions in terms of suppression (Raa <1) or enhancement (-Raa >1)- In contrast to 
the naive expectation the results for central Au+Au collisions at ^/snn = 200 GeV at 
RHIC show a suppression of up to a factor of five in the nuclear modification factor 
[3 Ej as shown in Fig. [TJ 
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Figure 2. Nuclear modification factor R&pS$>t) for neutral pions in d+Au collisions 
at -y/sNN = 200 GeV compared to theoretical predictions. The boxes around the data 
points indicate the systematic errors correlated in pt, and the shaded grey band at 1 
indicates the percent normalization error. 



It has been suggested that the observed suppression is a result of parton energy loss 
in a medium of high color-charge density (mainly gluons) like the quark-gluon plasma 
[2J E]. A number of model predictions (shown in Fig. ^) which were made prior to 
the release of the data can describe the data if they take into account constant (Wang 
et.al.) [Ij or energy- dependent (Levai et. al.) |Hj, (Vitev et. al.) [Oj parton energy-loss 
effects. Although the magnitude of the suppression in the data at pr ~ 4 GeV/c is 
qualitatively described by the model with constant energy loss, the increasing _Raa with 
Pt is contrary to what is observed in the data; whereas the energy-dependent energy-loss 
models give reasonable agreement with the data. 

In contrast to parton energy loss descriptions, a gluon saturation calculation ^U] 
is able to predict the magnitude of the observed suppression as well, though it fails 
to reproduce exactly the flat px dependence of the quenching. This model assumes a 
saturation of gluons at small Bjorken-x in the initial-state wavefunction of the Au nuclei, 
which leads to fewer hard gluon-gluon scatterings and thus to fewer high-pr particles. 
Similarly, estimates of final-state interactions in a dense hadronic medium ^1] yield the 
same amount of quenching as models based on partonic energy loss. So the latter two 
models claim to explain the observed suppression without assuming the formation of a 
quark-gluon plasma. 

Hence it is impossible to distinguish between initial-state effects (like gluon 
saturation) and final-sate effects (like parton energy loss) by studying solely Au+Au 
collisions at this energy. 
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Figure 3. Nuclear modification factor i?dA(pr) for neutral pions in d+Au collisions 
at ^/snn = 200 GeV for four different centrality selections compared to theoretical 
calculations. 



3. High-pr Particle Yields in gH-Au collisions at ^/snn = 200 GeV 

Explanations attributing the high-p^ hadron suppression in central Au+Au collision 
to final-state effects and models which are based on initial-state effects make different 
predictions for d+Au collisions: In a d+Au collision the highly excited matter is not 
created over a large volume and so high-p^ particle suppression is not expected in the 
parton energy-loss scenario. By contrast, in the framework of the gluon saturation 
model high-pr particle production ought to be suppressed in d+Au as well. The gluon 
saturation calculation in predicts -Rab ~ 0.7 for central d+Au collisions. d+Au 
reactions have been studied in early 2003 by PHENIX and first results have been 
published [TJ. 

The nuclear modification factor R^a for high-p^ 7r°'s in minimum-bias d+Au 
collisions is shown in Fig. El In contrast to central Au+Au collisions no suppression at 
mid-rapidity can be seen. Within uncertainties R^a for neutral pions is equal to unity 
in contrast to unidentified charged particles that exhibit a clear enhancement compared 
to p+p collisions (not shown). This enhanced production, observed first in p + A fixed- 
target experiments ( "Cronin effect" ) > is interpreted in terms of multiple initial-state 
soft and semi-hard interactions which broaden the transverse momentum of the colliding 
partons prior to the hard scattering. Figure El also shows that models incorporating both 
the Cronin enhancement and shadowing effects can describe the data well: A next-to- 
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Figure 4. Comparison of pion yields in the most central (0 — 20%) and the most 
peripheral (60 — 88%) data samples i?cp for three different measurements as a function 
of pt- 



leading-order pQCD prediction which was made prior to the release of the data (Levai 
et. al. [H]) with impact-parameter-dependent and independent shadowing effects and 
a model with a different implementation of Cronin and shadowing effects (Kopeliovich 
et. al. [13]) are compared to the data. These d+Au results support the conclusion that 
a final-state nuclear medium effect, such as parton energy loss, is necessary to describe 
the suppression observed in central Au+Au collisions. 

The influence and interplay of various nuclear effects on particle production depends 
on the geometry or centrality of the collision due to the changing nuclear density. 
Particle production at different rapidities also reflects the dynamics of the nuclear 
and Bjorken-x dependence of these effects. Therefore, a careful study of centrality- 
and pseudorapidity-dependence of particle production in d+Au collisions can provide 
an important baseline for the understanding of the relative contributions of different 
initial and final-state effects in Au+Au collisions. Here we will focus on the centrality 
dependence of neutral pion production - studies of pseudorapidity dependence can be 
found in [IB]. Figure El shows the nuclear modification factor RdA^Pr) for neutral pions in 
d+Au collisions at ^snn = 200 GeV for four different centrality selections compared to 
theoretical calculations. The boxes around the data points indicate the systematic errors 
correlated in p?, and the shaded grey band at 1 indicates the percent normalization error. 
Within systematic errors i?dA(pr) ~ 1- There is a hint for a small increase of i?dA with 
centrality. This would be in accordance with the observation for unidentified charged 
spectra which show a clear increase of RdA with centrality in d+Au, a trend opposite 
to the Au+Au results. 

The data can be well described over all centralities by a Glauber-Eikonal 
approach [TJ] shown as the red curve in Fig.Elwithout any assumption about initial-state 



6 




p T (GeV/c) 

Figure 5. Nuclear modification factor Raa{pt) for neutral pions in Au+Au collisions 
at ^/snn = 62.4 GeV for two different centrality selections compared to a theoretical 
prediction. 

gluon-saturation effects. The prediction for the gluon-saturation calculation (CGC, [TP] ) 
is plotted as well. Even with the systematic uncertainties of the preliminary PHENIX 
data the model prediction seems to underestimate the data in the most central data 
sample. 

To demonstrate the consistency of various measurements of charged and neutral 
pions within the PHENIX experiment Fig. 0] shows the ratio Rqp of pion yields in the 
most central (0 — 20%) to the most peripheral (60 — 88%) data samples each normalized 
by the respective (Tab) for three different measurements as a function of px- This 
ratio has the advantage that many of the systematic errors cancel. The figure shows 
measurements of charged pions at low px identified by the TOF detector [TH] and at 
high pt identified by the RICH and EMCal detectors together with the measurement 
of 7r°'s. All three independent measurements agree well within systematic errors. In 
future analyses it has to be proven if Rep for pt > 8 GeV stays at one. 

4. High-p T Particle Yields in Au+Au collisions at v /5nn = 62.4 GeV 

To study the onset of the suppression of neutral pions observed in central Au+Au 
collisions at ^/snn = 200 GeV and to map out the elastic and inelastic scattering 
properties of a possibly created quark-gluon plasma it is useful to study the a /snn 
dependence of neutral-pion production. To provide more information in the energy 
range between the measurements at CERN SPS ^Hj at a/snn = 17.3 GeV and the 
RHIC measurements at ^/s^ = 130 GeV and 200 GeV in early 2004 the RHIC 
experiments measured particle production in Au+Au collisions at -,/snn — 62.4 GeV. 
Figure El shows the results on the nuclear modification factor R\a{pt) for neutral pions 
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in Au+Au collisions at this energy for two different centrality selections. The boxes 
around the data points indicate the systematic errors correlated in p T , and the shaded 
grey band at 1 indicates the percent normalization error. Note the huge systematic 
error on the normalization due to large uncertainties in the p+p reference as there is 
no 7T° measurement from RHIC available at this energy [20] ■ Raa{pt) for 7r°'s shows 
a suppression of up to a factor 5 at high p T with a strong p T dependence. A model 
prediction which was made prior to the release of the data [21] incorporating the final- 
state energy loss with a gluon rapidity density of dN 9 /dy = 650 — 800 and initial-state 
Cronin scattering is shown in Fig. [5] as well. For an interpretation of the data a future 
measurement of the p+p reference at RHIC at this energy is strongly desirable. 



5. Conclusions 



One of the most significant observations at RHIC so far is the suppression of high- 
Pt hadrons at mid-rapidity in central a/snn = 200 GeV Au+Au collisions relative to 
binary-scaled p+p reference spectra. This observation is consistent with predictions in 
which hard-scattered partons lose energy in a dense medium. A competing theory, the 
gluon saturation model, appears as an unlikely explanation after the same measure- 
ments were made in d+Au collisions: Here no high-p^ particle suppression is observed 
at mid-rapidity supporting the expectation of the parton energy loss scenario but dis- 
proving predictions by the gluon saturation model. Very little centrality dependence of 
Raa(pt) for neutral pions is observed in d+Au. The i/snn dependence is studied by 7r° 
spectra in Au+Au collisions at ^/snn = 62.4 GeV. 

[1] PHENIX Collaboration, L. Aphecetche et at, Nucl. lustrum. Meth. A 499, 521 (2003). 

[2] PHENIX Collaboration, S.S. Adler et ai, Phys. Rev. Lett. 91, 241803 (2003). 

[3] PHENIX Collaboration, S.S. Adler et ai, Phys. Rev. Lett. 91, 072301 (2003). 

[4] PHENIX Collaboration, S. S. Adler et ai, Phys. Rev. C 69, 034910 (2004). 

[5] M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 (1990). 

[6] S.N. Wang, Nucl. Phys. A 715, 775 (2003). 

[7] X.N. Wang, Phys. Rev. C61, 064910 (2000). 

[8] P. Levai et ai, Nucl. Phys. A 698, 631 (2002). 

[9] M. Gyulassy, P. Levai and I. Vitev, Nucl. Phys. B 594, 371 (2001). 

[10] D. Kharzeev, E. Levin and L. McLerran, Phys. Lett. B 561, 93 (2003). 

[11] K. Gallmeister, C. Greiner and Z. Xu, Phys. Rev. C 67, 044905 (2003). 

[12] PHENIX Collaboration, S. S. Adler et ai, Phys. Rev. Lett. 91, 072303 (2003). 

[13] JW. Cronin et ai, Phys. Rev. D 11, 3105 (1975). 

[14] P. Levai, G. Papp, G. G. Barnafoldi and G. I. Fai, |arXiv:nucl-th/03060"T9| 

[15] B. Z. Kopeliovich, J. Nemchik, A. Schafer and A. V. Tarasov, Phys. Rev. Lett. 88, 232303 (2002). 

[arX"iv:hep-ph/020TuT0| . 

[16] PHENIX Collaboration, M. X. Liu, J. Phys. G 30, S1193 (2004). 

[17] A. Accardi and M. Gyulassy, Phys. Lett. B 586, 244 (2004). 

[18] PHENIX Collaboration, Adcox K et ai, Nucl. lustrum. Meth. A 499, 489 (2003). 

[19] WA98 Collaboration, M.M. Aggarwal et ai, Eur. Phys. J. C 23, 225 (2002). 

[20] D.d'Enterria, these proceedings. 
[211 



I. Vitev, arXiv:nucl-th/0404052 



